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ABSTRACT 

We analyze the cut-off dependence of mesonic spectral functions calculated at finite 
temperature on Euclidean lattices with finite temporal extent. In the infinite tem- 
perature limit we present analytic results for lattice spectral functions calculated 
with standard Wilson fermions as well as a truncated perfect action. We explicitly 
determine the influence of 'Wilson doublers' on the high momentum structure of the 
mesonic spectral functions and show that this cut-off effect is strongly suppressed 
when using an improved fermion action. 
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1 Introduction 



The thermal modification of basic properties of hadrons, e.g. their masses and 
decay widths, is one of the central issues in the discussion of experimental signals 
that can emerge from the dense partonic systems created in heavy ion collisions. 
Lattice calculations can, in principle, provide this information through the analysis 
of thermal properties of Euclidean correlation functions of suitably chosen operators 
carrying hadronic quantum numbers. These correlation functions contain all the 
necessary information on the temperature dependence of hadronic spectral functions. 
In order to make such studies quantitative and reliably extract information valid in 
the continuum limit we have to understand, however, in detail the cut-off dependence 
of spectral functions calculated on the lattice at finite temperature. We provide here 
a detailed analysis of mesonic spectral functions in the infinite temperature limit and 
discuss their cut-off dependence. This provides a basis for discussions of the cut-off 
dependence of spectral functions at finite temperature and is similar in spirit to 
studies of the cut-off dependence of the QCD equation of state which have first been 
performed in the ideal gas (infinite temperature) limit pQ. 

Information on the hadronic spectrum is extracted in lattice calculations from 
properties of Euclidean time correlation functions of suitably chosen hadronic cur- 
rents. It has been suggested recently [2] to apply the Maximum Entropy Method 
(MEM), a well known statistical tool for the analysis of noisy data [3], also to the 
analysis of these correlation functions. This opened the possibility to extract detailed 
information on hadronic spectral functions, a(u,T), at zero as well as finite tem- 
perature jUEl- The first studies of spectral functions based on the MEM approach 
[6-13] have indeed been encouraging. These first studies, however, also showed that 
it is necessary to get control over typical lattice problems like finite cut-off effects 
or the influence of fermion doublers on the energy dependence (u) of spectral func- 
tions. Cut-off effects show up in the large energy regime of spectral functions and 
one generally may not be too much worried about them when one is interested in 
extracting information about the low energy part of the spectral functions. This, 
however, is different for the analysis of properties of heavy quark bound states and 
also is of particular importance for studies performed at finite temperature. In the 
latter case the high energy part of, e.g. the vector spectral function, is directly re- 
lated to physically observable dilepton cross sections and has been studied in much 
detail in (resummed) perturbative calculations |14j . Moreover, in the plasma phase 
of QCD quasi-particle excitations are heavy and typically will have masses which 
rise proportional to the temperature. At high temperature this is expected to lead to 
broad resonance-like structures in spectral functions. In lattice calculations, which 
at present all have been performed with Wilson type fermion formulations, it is in 
general difficult to distinguish such effects from contributions arising from so-called 
heavy Wilson doublers. On a more technical level it also is important for the MEM 
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analysis to include information on the short distance behaviour of correlation func- 
tions in the default model. This requires information on lattice cut-off effects in the 
spectral function as well as possible modifications of the integration kernel. 

There are apparently plenty of reasons to get control over the lattice cut-off effects 
in numerical calculations of spectral functions. We will analyze these in the infinite 
temperature limit of QCD by explicitly calculating hadronic spectral functions on 
lattices with finite temporal extent N T . We will present results for spectral functions 
calculated on isotropic as well as anisotropic lattices and discuss their quark mass 
dependence. Moreover, we will present results for standard Wilson fermions as well 
as a truncated perfect action jT5J UH] • 

The paper is organized as follows. In Section 2 we summarize known results for 
free quark-antiquark spectral functions calculated in the continuum. In Section 3 we 
perform the corresponding calculation for Wilson fermions on the lattice. Section 4 
is devoted to a discussion of these lattice spectral functions, in particular their quark 
mass and anisotropy dependence. In Section 5 we present results from a calculation 
with an improved Wilson fermion action, a truncated perfect action. Finally we give 
our conclusions in Section 6. Some details of our calculations are presented in two 
Appendices. 

2 Thermal quark-antiquark spectral functions 

Thermal quark-antiquark correlation functions in coordinate space, Gjj(t,x), are 
defined as 



where (...) denotes the thermal average. The local sources for currents with different 
mesonic quantum numbers H are given by J#(t, x) = g(r, x)T H q(r, x), and T H is an 
appropriate combination of 7-matrices; i.e., Th = 1, 75, 7 M , 7^75 for scalar, pseudo- 
scalar, vector and pseudo-vector channels, respectively. From these we obtain the 
mixed correlation functions at fixed momentum p which are commonly considered 
in lattice calculations, 








These two point functions have a spectral representation, 




(2.3) 
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where an{uJ,p, T) denotes the temperature dependent spectral function and K(u, r) 
is the integration kernel which carries the entire dependence on Euclidean time 
re [0,1/T), 



cosh[w(r- 1/2T)] 
sinh(u;/2T) 



(2.4) 



It is easy to convince oneself that the spectral function appearing in Eq. 12.31 indeed 
is the Minkowski space spectral function (see e.g. [T7] ) 



-E n (p)/T, 



1-e 
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(2.5) 



where Z(T) is the partition function and (...) denotes here the matrix element of 
the hadronic current J#(0) taken between energy eigenstates at fixed momentum p. 

The correlation functions can be evaluated using the momentum space repre- 
sentation of the quark propagator and its spectral representation [T2j. In the fol- 
lowing we will mainly be concerned with the zero momentum spectral functions, 
(Jh(uj, T) = (Jh(u,p = 0, T), which take on a rather simple form [T7| ITS]. 
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Here iV c denotes the number of colours, eg. N c = 3. For some selected mesonic 
quantum number channels, H, the coefficients a#, bn and /# are given in Table l2~Tl 
We note that in some cases a contribution proportional to a 5-function at vanishing 
energy appears in Eq. 12.61 At non-zero temperature this gives rise to a constant, 
T-independent term in the Euclidean correlation function defined in Eq. 12.31 For 
massless quarks, also at non-vanishing momentum, p = a/p* 2 , a rather compact form 
for the spectral function is obtained, 
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Table 2.1: Coefficients a#, bjj and fu for free spectral functions in different mesonic 
quantum number channels H (Eq. 12. (jjl . The last two columns give coefficients 
appearing in the definition of the corresponding correlation functions on the lattice 
(Eq. 13. 5|) . The momentum dependent function d is defined in Eq. 13.71 

Due to asymptotic freedom of QCD the free field limit is approached at infinite 
temperature. In order to discuss the infinite temperature or free field limit of spectral 
functions and correlation functions it is appropriate to rescale all variables with non- 
trivial dimension by appropriate powers of the temperature, e.g. uj = uj/T . For fixed 
m = m/T the rescaled correlation functions, Gh{t, p) = Gh(tT, p/T)/T 3 , then have 
a well defined infinite temperature limit, 

POO 

G H (f,p)= dua H (u>,$,T) K(u,f) , (2.8) 
Jo 

where a denotes the rescaled spectral density, &h{&,P,T) = <jh(u,P,T)/T 2 , and 
k(Q,f) = cosh(w(f - l/2))/sinh(o)/2). 

Eventually we are interested in obtaining information on the temperature de- 
pendence of the physical spectrum and we thus want to determine the Minkowski 
space spectral function as introduced in Eq. 12.51 In the continuum the spectral func- 
tion is connected to the temporal Euclidean correlation function via Eq. 12.31 The 
analyticity properties required to make that connection, however, in general only 
exist in the continuum. Therefore, extracting a spectral function from lattice data 
on Euclidean correlation functions via Eq. 12.31 will suffer from lattice artifacts and 



5 



it is only in the continuum limit that a direct relation to the spectral properties at 
finite temperature can be established. This will become clear in the following section 
where we discuss the spectral representation of mesonic correlation functions using 
the standard Wilson fermion formulation on Euclidean lattices with finite temporal 
extent N T . 



3 Lattice Spectral Functions with Wilson Fermions 



In the following, the dimensionless representation of hadronic two point correlation 
functions and the corresponding rescaled spectral functions given in Eq. 12.81 will be 
analyzed for their cut-off dependence. The situation here is similar to the discus- 
sion of the cut-off dependence of bulk thermodynamic variables, e.g. the rescaled 
pressure, P/T 4 [Ud]. in that case, on isotropic lattices, deviations from the con- 
tinuum result can be expressed in terms of the lattice spacing given in units of the 
temperature, aT, which is nothing else but the inverse of the temporal extent of 
the lattice, aT = 1/N T . For bulk thermodynamic observables the temperature is in 
general the only a quantity with non-trivial dimension that can set the scale for the 
cut-off dependence. In the case of spectral functions, however, the energy provides 
another scale and we can expect to find an additional dependence of the spectral 
functions on au = uj/N t . However, we will show in the following explicitly that 
only Cj/N t determines the cut-off dependence and a sole dependence on 1/N T does 
not appear in the spectral functions. On anisotropic lattices cut-off effects are, in 
addition, controlled by the ratio £ = a/a T of spatial (a) and temporal (o T ) lattice 
spacings. In this case, temperature and energy in (spatial) lattice units also depend 
on this ratio as aT = £/N T and au = u£/N T , respectively. 



We will discuss here the cut-off dependence of spectral functions calculated 
within a generalization of Wilson's fermion discretization scheme [20] on anisotropic 
lattices. In the free field limit the fermion action is diagonal in the colour degrees 
of freedom, Sp = z Yll=i ^1 Qk,i with the fermion matrix 
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k— i.l 



(3.1) 



where the dimensionless quark mass, m, is expressed in units of the spatial lattice 
spacing, rh = ma. The generic choice of the Wilson action is r = r T = 1 20] • We 

a This is correct in the limit of vanishing as well as infinite quark masses. 
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will discuss here only the case r T = 1 as this avoids the occurrence of a second pole 
in the fermion propagator which would give rise to an unphysical time-like doubler 
mass. We will, however, consider the space-like Wilson r-parameter, r e (0,1], as 
an additional free parameter of the fermion action. Even on anisotropic lattices 
with £ > 1 the generic choice for the r-parameter is r = 1. Yet, by choosing 
r = l/£ discretization errors of order ma can be completely eliminated at leading 
order from meson correlators with improved Wilson (clover) fermions on anisotropic 
lattices |21j . e.g. in studies of heavy quark systems. As this choice has, in fact, been 
used in recent studies of heavy quark spectral functions [T3| we also will explore the 
dependence of spectral functions on the choice of r. 



Analytic results for the free field limit of hadronic correlation functions on 
isotropic lattices of size N% x iV T have been presented in the past using free Wil- 
son fermions [22| with r = 1. It is straightforward to extend these calculations 
to the case of anisotropic lattices and r / 1. Starting point for the calculation 
of hadronic correlation functions is the momentum space representation of the free 
Wilson fermion propagator, 



S(k) 
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(3.2) 
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(3.3) 
(3.4) 



On a finite lattice of size Njj. x iV r the momenta take on discrete values, ko = 
j^(n + 1/2) with n = 0, ± 1, .. ± (N T /2 - 1), N T /2 and h = j^rii with m = 
0,' ± 1, .. ± {N a /2 - 1), N a /2 for i = 1, 2, 3. 

Following Ref. 22 one finds for the temporal zero momentum free quark- ant iquark 
correlation functions with mesonic quantum numbers, H, 



G H (f,p = 0) = N c 



E 



cff(k) cosh[2E '(k)N T (f - 1/2)] + d l g(k) 
(1 + M) 2 cosh 2 {E{k)N T / 2) 



(3.5) 



which is defined on the discrete set of Euclidean times r = n/N T , with n = 
0, 1,...,N T — 1, accessible on a lattice with temporal extent iV T . The energy, 
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E = E(k), is given by the location of the pole of the denominator of the Wilson 
fermion propagator (Eq. 13.2)1 at ik = E(k), i.e. 



K? + M 2 

coshE = 1 + 2(TTM) • (3 ' 6) 

The functions c 1 ^ and d l jj appearing in Eq 13.51 depend on the three momentum k 
through the function 

"^Sb • (3 - 7) 

Note that d is approaching 1 in the continuum limit. For some quantum number 
channels explicit expressions are given in Table 12.11 

Eqs. 13. 5113. 71 can be used to analyze the infinite temperature limit of mesonic cor- 
relation functions on any finite lattice of size N% x N T . In the following we will take 
the thermodynamic limit b (N a — > oo) and concentrate on the cut-off dependence 
of these correlation functions which arises from N T being finite. In the thermody- 
namic limit the momenta are continuously distributed in the interval [— ir, 7r] and 
the energy consequently becomes a continuous function. This allows for an integral 
representation of Euclidean correlation functions in complete analogy to the contin- 
uum representation given in Eq. 12.31 In particular, we will show that also on lattices 
with finite temporal extent iV T the integration kernel is identical to the continuum 
kernel, Eq. 12.41 Cut-off effects which are responsible for the deviation of the lattice 
correlation functions from the corresponding continuum correlation functions thus 
only show up in the lattice spectral functions. 

In the thermodynamic limit the momentum sum appearing in Eq. 13.51 gets re- 
placed by a three-dimensional integral over the lattice Brillouin zone, 

k k -* 

It is obvious from Eq. 13.51 and Table 12.11 that there will appear only two types of 
r-dependent integrals which result from the two terms appearing in the definition 
of c#*, i.e. 1 and d{k) respectively, 



GAf) = N c — 



N T \ 3 f 1 cosh[2E(k)N T (f-l/2)] 



£ J J (l + M) 2 cosh 2 (E(k)N T / 2) 



3 In general we found that the dependence on the spatial extent becomes weak for £No-/N T >3. 
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A , ,N T \ 3 f d{k) cosh[2E(fc)iV r (f-l/2)] 



In addition there are also two r-independent integrals related to the sums involving 
the term d l jj appearing in Eq. 13.51 These constants contribute to the ^-functions 
at vanishing frequency. On finite lattices all quantum number channels listed in 
Table l2~T1 f except the pseudoscalar) will receive non- vanishing contributions f^, 

f iat = N (NA 3 [ df(k) 1 

lH C U J k(l + Mf C0 ^ { E{k)N T /2) ' 1 J 

The integrals given in Eq. 13 .81 are the starting point for our discussion of the cut- 
off dependence of mesonic spectral functions. So far we have achieved to represent 
the correlation functions Gh(t) in terms of three-dimensional momentum integrals. 
Our goal is to find, for finite iV r , a spectral representation defined through the 
one-dimensional integral given in Eq. 12.31 This can be achieved by introducing the 
energy in units of the temperature, Q = 2EN T , as one of the integration variables. 
In Appendix A we show explicitly for the case r = 1 the sequence of variable 
transformations required to obtain an integral representation which is in complete 
analogy to the continuum relation, Eq. 12.81 i.e. we can write the integrand of this 
integral as a product of a r-independent lattice spectral function expressed in units 
of T 2 and a r-dependent kernel which is identical with the continuum kernel K 
defined in Eq. EU 



Gi(f)= / dCo d^Co^Nr) K(u,f) , i = 1, 2 . (3.10) 



The spectral functions, N T ), explicitly depend on the lattice cut-off, which is 
reflected in the explicit dependence on u/N T = ua T , 



UO \ I UJ \ I UJ 



N T ) = tanh (J j sinh j sinh ) h(uj/N T , £) 

at at2 / ~.\ sinh 3 ( jfir 
a 2 (Co,N T ) = tanh - I a (u>/N T ,£) . (3.11) 

The two-dimensional integrals Ii(u/N T ,£) are worked out in Appendix A for the 
case r = 1. They can, however, also be defined for arbitrary values of r. 
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The integration limits in Eq. 13.101 depend on the quark mass, the anisotropy and 
the Wilson r-parameter. For r = 1 the maximal energy is determined by the largest 
quark three momentum possible, k = (it, tt, 7r), and we find from Eq. 13. fit 

u min = 2N T ln(l + m/0 , u) max = 2iV T ln(l + (6 +m)/f) . (3.12) 

For r < 1, however, the maximal energy generally corresponds to a momentum in 
the interior of the first Brillouin zone. The corners of the 3-dimensional Brillouin 
zone are local minima of the dispersion relation which are interpreted as doubler 
masses which are proportional to the Wilson r-parameter, i. e. they become lighter 
with decreasing r. As we will see this leads to rather complicated spectral properties 
even in the free quark, infinite temperature limit. 

We also note that the lattice spectral functions cr, are directly proportional to the 
massless spectral functions in the continuum (&h ~ a) 2 tanh(tu/4)). In the massless 
limit the deviations from the continuum results thus only arise through the ratio 
Cj/N t = ua T which is the energy expressed in units of the temporal lattice spacing. 
Eq . 13 . 1 1 1 explicitly reflects a well known feature of the lattice formulation, i.e. cut-off 
effects depend on the energy scale. Of course, an explicit dependence on the lattice 
spacings (a, a T ) never appears in the lattice formulation. It is, however, remarkable 
that also no explicit dependence on 1/N T = a T T appears in the spectral functions. 

Finally, we want to reconstruct from <7j the spectral functions in fixed quantum 
number channels. In particular, we will consider spectral functions in the pseudo- 
scalar, scalar, vector and axial-vector channels which are given by 

a l p S = <7i , ag* — — cr 2 , b l y = 3<7i - a 2 , a a — a\ — 3<t 2 . (3.13) 

We ignore here a term proportional to Cj8(uj) which, as discussed above, arises from 
the r-independent part (/^*) of the correlation functions. For the above quantum 
number channels the coefficients of the (^-functions will vanish in the continuum limit 
and we have checked that they are indeed small already on lattices with temporal 
extent N T ~ 10. 

4 Cut-off effects on isotropic and anisotropic lat- 
tices 

We will analyze here in detail the lattice spectral functions derived in the previous 
section for some choices of paramters (r, £, m) which have been used in recent 
studies of meson spectral functions at finite temperature. 
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4.1 Massless Quarks on isotropic lattices: £ = 1, m = 



Let us first discuss the lattice size dependence of the spectral functions for the 
case of vanishing quark masses (m = 0) and on isotropic lattices (£ = 1). In 
Fig. 14. If left) we show the ratios of the lattice and the corresponding continuum 
spectral functions. We note that this ratio is a function of u/N T = ua only. The 
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Figure 4.1: Ratio of lattice and continuum spectral functions with m = and £ = 1 
and for different quantum number channels (left). The figure on the right hand 
side shows the pseudo-scalar lattice spectral functions calculated for lattices with 
temporal extent N T = 16, 24 and 32. The solid line corresponds to the continuum 
spectral function given in Eq. 12.61 



spectral functions vanish for u > tu max , where u mSuX /N T = u max a = 2 In 7 is obtained 
from Eq. 13.61 as the largest energy possible for a mesonic state constructed from 
two independent massless Wilson fermions with momentum k = (tt,tt,tt). In the 
pseudo-scalar and vector channels we observe a pronounced peak which occurs when 
the momenta of both Wilson fermions correspond to the first corner of the Brillouin 
zone k = (tt, 0, 0). The corresponding energy is u%a = 2 In 3. Finally we observe a 
cusp at u-ia = 2 In 5 between u±a and u metK a which corresponds to the second corner 
of the Brillouin zone k = (rr, tt, 0). In the continuum limit the lattice artifacts shift 
to higher energies and the lattice spectral functions approach the continuum result 
(Fig. EIKright) ). c 

Also in the interacting case [Uj a peak-like structure has been observed at similar 
values of the energy, i.e. for ua ~ 1.7 in the pseudo-scalar and ua ~ 2 in the 
vector spectral function, respectively. It has been shown that these peaks shift to 
larger energies when decreasing the lattice spacing. They thus have been identified 
as lattice artifacts which do not correspond to physical states in the continuum limit 
jH]. It is likely that these peaks are remnants of the cut-off effects seen here in the 

c Note that in the scalar and axial case we always plot the positive functions — as and — a a- 
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free spectral functions. To establish this relation in more detail it would certainly be 
interesting to analyze also spectral functions in the scalar and axial- vector channels. 
At least in the free case the peaks which show up in the pseudo-scalar and vector 
spectral functions are absent in the scalar and axial-vector channels. 

In Ref. 6 it has been suggested that the cut-off dependent peaks in the spectral 
functions are related to bound states involving heavy quark doublers with masses of 
0(1/ a). In view of the free spectral function we would, however, prefer not to speak 
of states at all. Rather, the distortion of the spectral function and the characteristic 
structures seen in the free case are due to the lattice dispersion relation and to the 
sudden restriction of available momentum space that occurs when one of the fermion 
momenta reaches one of the corners of the first Brillouin zone. 

We also note that all four spectral functions coincide up to wa ~ 1.5 where 
they differ by less than 15% from the continuum result. This agreement of different 
quantum number channels is reminiscent of the chiral symmetry of the free fermion 
action. In the lattice formulation with Wilson fermions chiral symmetry is, of course, 
explicitly broken which leads to different spectral functions in the scalar and pseudo- 
scalar sector. As can be seen from Fig. 14.11 this explicit breaking most strongly 
influences the large energy part of the spectrum, i.e. aps deviates strongly from 
— eg for ua > 1.5. The same holds true in the vector and axial- vector channels. Also 
note that for large energies the finite cut-off effects can lead to a negative lattice 
spectral function in the axial-vector channel. 



4.2 Massless Quarks on anisotropic lattices: £ > 1, m = 

When reducing the temporal relative to the spatial lattice spacing (£ = a/a T > 1) 
one has to increase the number of grid points in the temporal direction if one wants to 
keep the temperature constant, 1/T = N T a T . In the interacting case it requires a fine 
tuning of spatial and temporal couplings (hopping parameter) in order to maintain 
rotational symmetry at zero temperature and, of course, it will also increase the 
computational effort. Nonetheless, it may be of advantage in the analysis of mesonic 
correlation functions at high temperature because one can make use of information 
on the correlation functions at a larger number of Euclidean time steps. 

r = 1 : 

Let us first discuss the structure of free spectral functions on anisotropic lattices 
for the case r = 1. As can be seen from Eq. 13.121 for a fixed ratio £/N T = Ta the 
support for the spectral function increases with increasing anisotropy. For (N T , £) — ► 
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oo it reaches a finite limit, i.e. uj m i n — > 2fh and Cj max —>■ 12N T /^ + 2m. By increasing 
the temporal lattice size and the anisotropy simultaneously the upper limit, above 
which the spectral functions vanish, can thus be increased by about a factor 3 relative 
to the case of isotropic lattices. For a moderate anisotropy factor (£ ~ 4), typically 
used in numerical calculations, £b max (and thus also oo max a ) is about a factor 2 larger 
than in the isotropic case. Lattice artifacts, however, set in earlier; for £ = 4 the 
peak in the pseudo-scalar and vector spectral function is only shifted by a factor 1.4 
(see Fig.lOJ). 

In Fig. l4.2l we show the same spectral functions as in Fig. l4.1l now calculated with 
an anisotropy £ = 4; a choice of the anisotropy parameter which has been used in 
recent studies of spectral functions ^21 IS] ■ As can be seen from this figure the en- 
ergy interval in which the pseudo-scalar and vector spectral functions are only little 
affected by cut-off effects is about twice as large as on the isotropic lattices. Scalar 
and pseudo-scalar correlation functions, however, are affected differently. As a con- 
sequence the degeneracy of both spectral functions, a precursor of chiral symmetry 
restoration, is lifted at smaller energies than on isotropic lattices. The situation is 
similar for the vector and axial- vector spectral functions. 



2.50 
2.00 
1.50 
1.00 



0.50 



0.00 



£=4 /\ 


v 

PS - - 
s — - 

A 









0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 



Figure 4.2: Ratio of lattice and continuum spectral functions for m = and £ = 4. 



r = l/£: 



When one uses the particular choice r = l/£ jTS| the doubler masses become 
light for large £ and can influence the spectral properties at lower energies than it 
is the case for r = 1. As the quark dispersion relation no longer leads to maxima 
in the corner of the Brillouin zone the energy range in which the spectral functions 
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Figure 4.3: Pseudo-scalar lattice spectral function for massless quarks calculated 
on a lattice with temporal extent N T = 96 on an anisotropic lattice with £ = 4 
and Wilson r-parameter r= l/£. The dashed curve shows the sum of contributions 
from continuum spectral functions for one massless pseudo-scalar and a number of 
massive pseudo-scalars constructed from free quarks with masses corresponding to 
the doubler masses obtained from the quark dispersion relation of this action. The 
horizontal lines indicate the location of doubler masses in the three distinct corners 
of the first Brillouin zone. 



are non-zero shrinks compared to the r = 1 case. For £ = 1/r = 4 one finds from 
the quark dispersion relation uj max a = 3.45, which is even slightly smaller than the 
corresponding value for r = £ = 1. 

The non-monotonic behaviour of the dispersion relation also makes it more com- 
plicated to find a closed analytic representation for the spectral functions in terms of 
two-dimensional integrals as we have done for the case r = 1. Although in principle 
it is possible to generalize the approach described in Appendix A for r = 1, in the 
case r < 1 we have used the simpler numerical binning-approach to calculate spec- 
tral functions. This is also introduced in Appendix A. The resulting pseudo-scalar 
spectral function is shown in Fig. 14.31 for massless quarks and £ = 1/r = 4. In this 
case the mass of the lightest, three-fold degenerate doubler is m\a = 0.471 which 
gives rise to the first threshold at u\ = 2N T m\a/^ seen in this figure. The other 
structures seen in this figure result from contributions of the doublers in the other 
corners of the Brillouin zone as well as the maxima of the quark dispersion relation 
which now do not reside in the corners of the Brillouin zone. 

We note that this action does reproduce the continuum spectral function well up 
to the point where the first doubler starts contributing, uoa ~ 1. The energy range 
in which a good agreement with the continuum spectral function can be achieved 
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thus is compatible with the isotropic case. 



4.3 Massive Quarks 

The modification of heavy quark bound states and, in particular, their dissolution in 
a quark gluon plasma is considered to be an important signature for the formation 
of dense matter in heavy ion collisions. It thus is of interest to analyze also spectral 
functions for heavy quark bound states in lattice calculations. The first attempts to 
do so followed the strategies discussed in the previous sections, i.e. calculations have 
been performed with Wilson or clover fermions on isotropic ^T] and anisotropic ^3] 
lattices. In the latter case the formulation with r = l/£ has been used. 

The analysis presented in the previous sections for massless quarks carries over 
also to the case of massive quarks. The structure of the cut-off dependence discussed 
for the various types of actions follows patterns similar to those seen in the massless 
case. In the large energy region the cut-off effects dominate and a non-vanishing 
quark mass has little influence on the location of the pronounced peaks observed 
for r = 1 or the additional thresholds arising for r — l/£. This is shown in Fig. 14.41 
where we present results for quark masses m/T = 4.8. Aside from the low energy 
threshold which now is shifted to uj ~ 2m /T the spectral functions are similar to 
those shown in figures 14.11 and 14.31 

5 Lattice Spectral Functions with improved Wil- 
son Fermions 

In the previous sections we have seen that spectral functions obtained from calcula- 
tions with the standard Wilson action reproduce the continuum spectral functions at 
low energies, u;a<1.5. Deviations from the continuum dispersion relation, however, 
lead to strong modifications of the spectral functions at larger energies. Moreover, 
the violation of chiral symmetry becomes visible in spectral functions at these en- 
ergies. This motivates to analyze the cut-off effects of hadron correlation functions 
computed with improved fermion actions which have better chiral properties and 
lead to an improved dispersion relation. Much progress has been made in construct- 
ing such actions j2H|. As a first step in this direction we want to analyze here a 
simple truncated version of a fixed point action f^lCni- This action is constructed 
from a small set of couplings which connect sites in an elementary hypercube of the 
lattice and can be handled in close analogy to the case of the Wilson action. We 
restrict ourselves to the discussion of the massless case on isotropic lattices. The 
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Figure 4.4: Free pseudo-scalar spectral functions for fixed m/T calculated on 
isotropic lattices with different temporal extent and quark masses. Also shown 
is the result obtained on an anisotropic lattice with £ = 1/r = 4. In this case the 
spectral function deviates strongly from the continuum result for u)>30 (see also 
Fig l4.3|) . For better visibility we have cut-out this part here and replaced it by a 
wavy line. 
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action can be written as 



s = i 'vv^ -y) + K x - y) mi/) 

x,y V=0 



(5.1) 



with 



pA x - y) 



+ Yl P3(^y,x+fi+u+p — <5y,x~fi+u+p) 
+ ^ ] Pi(&y,x+p+D+p+cr ~~ $y,x 



',x— fi+v+p+a ! i 



(5.2) 



X(x - y) 



~T~ "y,x—[i+v j 
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',x—fi+v+p) 



+ ^(5y,x+fi+u+p+& + Sy,x-[i+D+p+6) ■ (5.3) 



Here fi, z>, p and <r denote unit vectors along positive as well as, except for /t, negative 
directions in the hypercubic lattice. Numerical values for the set of nine couplings 
{pi}f =1 , {Xi}t =0 are given in Table 1 of [TE] for m = 0. 

Taking the Fourier transform of the action, Eq. (J5.1)) . it is straightforward to 
write down the propagator in momentum space 

, . _ -«7o<5 sin ko — i^i ~ iK>i cos k + k x + k 2 cos k 

{ ' ~ (/C 2 + k\ + 5 2 ) + 2 cos fc (/Ci/C 2 + Ki k 2 ) + cos 2 k (JC 2 2 + k 2 2 - 5 2 ) ' 



with 



3 3 
/Ci = ^7i«i , K, 2 = y^7iA • ( 5 - 5 ) 



i=l »=1 



Explicit expressions for the momentum dependent functions ctj = a* (A;), $ = 
Pi(k), 5 = 6(k), K\ = Ki(k) and k 2 = k 2 (A:) are given in Appendix B. 

For the analysis of meson correlation functions and their spectral representation 
one has to calculate the quark propagator in the mixed (r, k) representation. The 
calculational steps are completely analogous to those for the standard Wilson action 
presented in Ref. 22. However, the quark propagator now has two poles, ik = E{. 
This is similar to the case of the Wilson action with r T < 1. The two poles are 
determined from 

coshEl = -_^wwEm, ( 5.6) 

(0\ Vi c 1 -P + y/(P 2 - ggj 
sgn(Q) cosh E 2 = — , (5.7) 



with additional functions 



P(Jfe) = /C1/C2 + «i«2, Q(k) = K% + k\ - 5 2 , R{k) = K 2 1 + K 2 l + 5 2 . (5. 
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Figure 5.1: The physical (E±) and unphysical (E%) poles of the quark propagator 
obtained from the truncated perfect action as the function of the spatial momentum 
k. The straight line shows the continuum dispersion relation and the dash-dotted 
curve gives the dispersion relation for the standard Wilson fermion action (E w ). 



In the limit N T — > oo the quark propagator is then given by 

Soo(r,k) = — 1 (5.9) 

[(«], — z/Ci) + (k 2 — iK, 2 ) cosh Ei + 74 5 sgn(r) sinh E^] e - ' 5 ' 17 " 



2^P 2 -QR sinh£ 2 

[(«i - z/Ci) sgn(Q) + («2 - ^2) cosh £ 2 + 74 <*> sgn(r) sinh £ 2 ] e" B2T . 

The first term describes the propagation of a physical state with the energy (disper- 
sion relation) Ei(k) while the second corresponds to an unphysical state, the analog 
of the time doubler in the case of the Wilson action with ry < 1. In Fig. 15.11 we 
show Ei and E 2 as function of k. 

As one can see from the figure E\ (k) is very close to the continuum result for small 
and moderate momenta and E 2 is much larger than Ex. Only for momenta \k\ > 2.5 
the gap between E 2 and Ei is getting small. This has important consequences for 
the meson correlators which will be discussed below. Following exactly the same 
procedure as for the Wilson action the meson correlators for finite N T can be written 
as 

G (fS=0) = N ( Nt ) 3 ^ °" ^ cosh i 2E ^)N T (f - 1/2)] + dff*(fc) 
c[ Nj Y (P 2 -QR)cosh 2 (Ei(k)N T /2) 

+ AGj(f). (5.10) 
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dP 




35 2 - dP 


dP-5 2 




5 2 -d? 


dP 


i=i 


-2d? 


2d p - 35 2 


M=0 


5 2 - 3d p 


3(d p - 6 2 ) 



Table 5.1: The coefficient c l ^ t,p and d l ^' p appearing in Eq. 15.101 The functions d p 
and S are defined in Eq. 15.111 and Eq. IB .71 respectively. 

The functions c?f' p {k) and d l ^ t ' p (k) are given in Table IB~T1 in terms of 5(k) and a new 
function d p (k), 



d p = d p {k) 



(Xi + JC 2 cosh E 1 
sinh 2 Ei 



(5.11) 



The term AGf (r) in Eq. 15.101 contains contribution from the second pole E 2 and 
an explicit expression for it is given in Appendix B. As the energies E 2 are large it 
turns out that this part leads to negligible contributions to the correlation functions 
except for very short distances (f = 0, 1/N T ). Using the binning-approach discussed 
in Appendix A we have calculated numerically the spectral functions in different 
quantum number channels. Results for N T = 24 are shown in Fig. 15.21 The good 
chiral properties and agreement with the continuum result over a wide range of 
energies (ua ~ 5) is self-evident. The contribution of the physical pole (E p ) has also 
been calculated analytically, which can be done in complete analogy to the case of 
the Wilson action. This contribution alone is, as expected, undistinguishable from 
the complete result up to ua ~ 5, Fig. 15.21 



19 



0.20 


~ lat/~\ /- 2 
■<7 H (W)/W 


V/2 

i'i PS - — 


0.15 
0.10 


iV r =24 


L\ s 

! \\ A/2 

j '.'font 

i ;', 


0.05 






0.00 


(_ 


! \.„_ w 




— I — — 


0.05 







- lat/~N ;~ 2 
o-jj (w)/w 



iV=24 



V/2 
PS 

s 

A/2 
cont 



\ '■'*■• 



20 40 60 80 100 120 140 160 180 



20 40 60 80 100 120 140 160 180 



Figure 5.2: Free spectral functions on isotropic lattices calculated with a truncated 
perfect action for m — 0. The left figure shows the complete result, the right one 
the contribution from the physical pole (Ef) alone. 



6 Conclusions 

We have presented an explicit calculation of mesonic spectral functions in the in- 
finite temperature limit of lattice QCD. We have analyzed the cut-off dependence 
of spectral functions in different quantum number channels for the Wilson fermion 
action on isotropic as well as anisotropic lattices. We find that the cut-off effects 
are of similar magnitude in both cases. The introduction of a Wilson r-parameter 
being less than unity in particular the choice r = l/£ does not seem to lead to a 
significant reduction of cut-off effects. 

Furthermore, we have analyzed the spectral representation of mesonic correlation 
functions using a truncated perfect action. As expected this does lead to a drastic 
improvement of spectral functions; cut-off effects are shifted to the high energy 
regime and chiral symmetry is preserved in the low energy part of the spectral 
functions up to energies ua ~ 5. 
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A Appendix 



We derive here the explicit form of the representation of the free mesonic lattice 
spectral functions given in Eq. 13.111 In particular, we give closed analytic expression 
for the integrals Ii(y,£) appearing in these equations for the case r = 1 and discuss a 
representation in terms of finite Riemann sums which is more convenient for dealing 
with the case r < 1 or spectral functions resulting from more complicated actions 
like the truncated perfect action analyzed in Section 5. 



1: 



We start from Eq. 13.81 by performing the obvious variable transformation k — > 
x= {sm 2 {k 1 /2),sin 2 {k 2 /2),sm 2 {k 3 /2)). We then define 

IC 2 + M 2 



and rewrite Eq. 13.61 as E = 2\sx{\/a + y/\ + a). Using these expressions we can 
perform a second variable transformation, x — > (uj,z 2 ,z 3 ) = (2N T E, x 2 /ot, x 3 /a). 
This leads to the representation of the correlation functions Gi (Eq. 13.10(1 where the 
integrals Ii(y,£) (Eq. 13.11]) are given by 



li(y,0= dz 2 dz 3 A(z) B(z) Q(z) , (A.2) 

JjQ(y) 

with z — (z,z 2 ,z 3 ) and z(y) = sinh 2 (?//4). We further define 
e - (Z2 + z 3 )(l - 2£z + m) + £m - f - 2zz 2 z 3 

X = ; : • A. 3 

1 + 2z(z 2 + z 3 )+m-2£z v ; 
With this the functions appearing in the integrand of Jj can be written as 

A(z) = e ^ + 2 ^ + ^ + ^3)+m 

y J ^ 1 + 2z{z 2 + z 3 ) + rh - 2& ' V ' ; 

B(z) = (l + 2C 1 z{x + z 2 + z 3 )+C 1 my 2 , (A.5) 

d(z) = (x z 2 z 3 (1 - zx) (1 - zz 2 ) (1 - zz 3 ))- 1/2 , (A.6) 
r] 4 x(l - zx) + z 2 (l - zz 2 ) + z 3 (l-zz 3 ) 

2[Z) " e{xz 2 z 3 {l-zx){l-zz 2 ){l-zz 3 ))^ ' { ■ } 
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The boundary of the two-dimensional integration region, fl(y), is given by 

n(y) = {z 2 , z 3 \0<xy<l; < z 2 y < 1; < z 3 y < 1} . (A.8) 

arbitrary r: 



In general we can determine the spectral function in a given quantum number 
channel also directly from the representation of the correlation functions in Eq. 13.51 
We can divide the interval of non- vanishing energies, Co = 2N T E(k), in n bins. 
Denoting by cjj the central values of these bins and introducing the bin length, 
e = {Cj max — tj min )/n, we can approximate Eq. 13.51 bv 



n 



cosh(£;(f- 1/2)) 



G H (f,p = 0)~^ e Si : " . ' " , (A.9) 

^ smh(o; i /2) 

where s« receives contributions from all terms in the momentum sum which lead to 
energies 2N T E(k) in the i-th bin, 



N c (J^j J2 ®( 2N rE(k) - (n - l)e) &(ne - 2N T E(k)) ■ 

c l g(k) smh(N T E{k)) 



[1 + M) 2 cos\i 2 {E(k)N T /2) 



(A.10) 



In the limit of large spatial volumes (N a — > oo) and small bin sizes (e — * 0) this 
gives the spectral functions in a given quantum number channel. We have used this 
approach to calculate spectral functions in the case r < 1 and also in the case of the 
truncated perfect action. Typically we used N a ~ 1000 and n = 1000. 



B Truncated perfect action 

In this Appendix we give explicit expressions for the auxiliary functions which ap- 
pear in calculations with the truncated perfect action, in particular in the quark 
propagator given in Eq. 15.41 
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The explicit form of the functions ai(k), (3i(k), % = 1,2,3 and 5(k), Ki(k) and 
K 2 (fc) are given below 

ai(fe) = 2si0oi + 2p 2 (c 2 + cs) + 4p 3 c 2 c 3 ) (B.l) 

a 2 (fc) = 2s 2 (p! + 2p 2 (c! + c 3 ) + 4p 3 c 1 c 3 ) (B.2) 

a 3 (fc) = 2s 3 (pi + 2p 2 (ci + c 2 ) + 4p 3 cic 2 ) (B.3) 

= 4sx(p 2 + 2p 3 (c 2 + c 3 ) + 4p 4 c 2 c 3 ) (B.4) 

(3 2 {k) = 4s 2 (p 2 + 2p 3 (c! + c 3 ) + 4p 4 C!C 3 ) (B.5) 

fo(k) = 4s 3 (p 2 + 2p 3 (ci + c 2 ) + 4p 4 cic 2 ) (B.6) 

S(k) = 2pi + 4p 2 (ci + c 2 + c 3 ) + 8p 3 (cic 2 + c 2 c 3 + cic 3 ) + 16p 4 cic 2 c 3 (B.7) 

rei(A;) = A + 2Ai(ci + c 2 + c 3 ) + 4A 2 (cic 2 + c 2 c 3 + cic 3 ) + 8A 3 cic 2 c 3 (B.8) 

« 2 (fc) = 2Ai + 4A 2 (ci + c 2 + c 3 ) + 8A 3 (cic 2 + c 2 c 3 + C1C3) + 16A 4 cic 2 c 3 , (B.9) 



where we introduced the shorthand notation q = cosfcj and = sinfcj. 

Furthermore, we give the explicit expression for the term in the meson correlator 
(f) coming from contributions of the second pole E 2 . It can be written as 

AGh 2 (t) — G H2 {f) + G m2 {f) , (B.10) 

where Gh 2 (t ) contains contributions from the second pole ik = E 2 only and there- 
fore can be written down in close analogy with the term containing only the contri- 
bution from the first pole E%, 

r ifa-m v^M'V c«h[2g 2 (S)iV r (f - 1/2)] + d'g-'Hk) 
G H2 (r,p = 0) = iV c (-) ^ ———^—^ .(B.ll, 



The functions P, Q, R have been defined in Eq. 15.81 d l ^ t,p2 and c l ^ t,p2 have exactly 
the same structure as d^' p and c^' p listed in Table l5~T1 with d p being replaced by 

^^^(fe) = (/Cl + /C2Sgn l g)CQShE2)2 , (B.12) 

sinh E 2 

The second term G//i 2 (r) contains contributions from both the first and second pole 
and can be written as 
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2df - Adf 
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Table B.l: The explicit form of the functions g^' p appearing in Eq. IB. 131 The 
functions df(k) and df(k) are defined in Eqs. IB. 14l and IB~15l 



GWf) = N c 



N a J ^ (P 2 - QR) cosh(Pi(£)iV T /2) cosh(P 2 (£)iV r /2) 



g l ^ p (k){cosh[E s (k)N T (f - 1/2)] - cosh[E d (k)N T (f - 1/2)]} 
+ 8 2 (k){cosh[E s (k)N T (f - 1/2)] + cosh[E d (k)N T (f - 1/2)]} 



(B.13) 



where E s = E\ + P 2 , E d = E\ — P 2 . An additional function g l ^ t,p {k) has been 
introduced which depends on 

^ P 3 _ fc-i sg n (<5) + cos h ^1 cosh E 2 + /Ci/C 2 (cosh P 2 + sgn(Q) cosh Pi) ^ 
1 ~ sinh E x sinh P 2 ' 



^ p3 _ fi^ sgn(Q) + k\ cosh Pi cosh P 2 + KiK 2 (cosh P 2 + sgn(Q) cosh E x ) 

sinh Pi sinh P 2 



The explicit form of g l ^ t,p for different quantum numbers is given in Table IB.l 
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